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external harmonic confinement. Using an efficient Green's function method we evaluate 
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gas with point-like hard-core interactions. 
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1. Introduction 

Atomic Fermi vapours have been cooled down to degeneracy in several experiments on 
ultra-cold atoms of 6 Li and 40 K [ffl 0, ^, using techniques similar to those employed 
to obtain Bose-Einstein condensation M. At very low temperatures collisions between 
atoms are essentially in the s-wave channel; however, all even-channel collisions are 
forbidden for spin-polarized fermions by the Pauli exclusion principle. Since dipolar 
interactions are also negligible, the Fermi gas can be considered as non-interacting. In 
the experiments the atoms are trapped by magnetic or optical potentials, providing a 
confinement which can be well approximated as harmonic. The geometry of the trapping 
potential can be very anisotropic, thus allowing to reach quasi-lD or quasi-2D situations 
i- 

We focus here on the quasi- ID geometry, which is the case of the Paris experiment 
and is relevant for experiments on atom lasers and atomic waveguides. The ID 
gas of non-interacting fermions serves also as a model for describing the ID Bose gas 
with hard-core point-like interactions (Tonks gas, Indeed, it was demonstrated by 
Girardeau that there exists an exact mapping between the bosonic and the fermionc 
many-body wave functions, = \^f\, valid also in the time domain |7|, ||. The Tonks 
gas limit corresponds to the case of total reflection in a two-atom collision ||, and the 
interactions in ID mimic the role of the Fermi statistics. 

In this paper we present some results for the exact equilibrium profiles and for the 
dynamical properties of a ID Fermi gas in harmonic confinement. In Sec. |2| we illustrate 
the role of quantum statistics on the particle and kinetic-energy density profiles, and we 
investigate how the Pauli principle influences the pair correlation function. In Sec. |3] we 
evaluate the spectrum of the gas both in the collisionless and in the collisional regime. 
In addition, we relate the result for the frequency of the quadrupole mode to a scaling 
property of the Hamiltonian used in the long-wavelength limit. Finally, Sec. [| offers a 
summary and some concluding remarks. 

2. Equilibrium properties 

We consider a system of N non-interacting Fermi particles confined in an external 
potential which is very tight in two spatial directions. We assume that the gas is in 
a quasi- ID situation where the transverse degrees of freedom are frozen, that is, the 
chemical potential is much smaller than the level spacing in the transverse direction. In 
this case the transverse wave function is the ground state one, which will factorize out 
from all the following calculations. We can thus focus only on the axial coordinate x 
and define the external confinement along x given by a real potential V ex t(x). For the 
sake of simplicity we will assume that the eigenfunctions ipi(x) are real ||10|| . 

The equilibrium properties of the gas can be expressed in terms of the longitudinal 
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one-body density matrix, which at zero temperature reads 

N 

^(x,z') = X)^WiM- (1) 

1=1 

In the following we will give several examples of calculations of the exact spatial profiles 
obtained as the moments of the density matrix, and we will evaluate also the longitudinal 
pair distribution function. 

2.1. Density profiles 

The zeroth moment of the one-body density matrix yields the equilibrium density profile 
n(x) 

n{x) = p^(x,x')\ x > =x = \ S i x ~ x i)\i>i), ( 2 ) 

i 

where are the positions of the fermions and | ipj) are the eigenstates of the Hamiltonian 
in the coordinate representation. We have shown [TT| that n(x) can be written in terms 
of the Green's function in coordinate space, G(x) = (x — x + ie)~ l , which is a function 
of the position operator x: 

1 N 

n{x) = lim Im / vV'i 1^(^)1 V'i) ■ (3) 

7T £^0+ L — ' 
i=l 

This equation expresses the density profile as the trace of the operator G(x) on the first 
N quantum states (Tr/v). Using the relation between the partial trace Tr^r of a generic 
matrix Q and the determinant of the inverse matrix Q -1 , 

Tr N Q = 9[lndet(g- 1 + XI N )]/d\\ x=0 , (4) 

we obtain 

1 d 

n(x) — lim Im — -[lndetfx + is — x — AIat)1a=o • (5) 

it £^o+ d\ 

This alternative expression for the density profile is very practical for performing 
numerical calculations in the case of harmonic confinement, V ext (x) = mujf w x 2 /2. For 
this specific potential the representation for both the position and the momentum 
operators on the basis ipi of the eigenstates of the harmonic oscillator is a tridiagonal 
matrix. Thus the operator x + ie — XIn is tridiagonal: the determinant of such a matrix 



can be evaluated recursively through a renormalization procedure pT| , |12| . 

The numerical results for the particle density profiles are shown in Fig. [1] at 
various number of fermions, as compared with the predictions of the local-density 
approximation: the latter neglects the shell structure and the spill-out from the classical 
border. Because of the symmetry between position x and momentum p in the harmonic- 
oscillator Hamiltonian, the momentum distribution has exactly the same behaviour as 
the particle density illustrated in Fig. |I[ 

On account of the mapping with the Tonks gas the particle density profile of iV 
noninteracting fermions in the ID harmonic trap is identical to that of iV bosons with 
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hard-core interactions in the same external potential, while the momentum distribution 



for such a system has a sharp peak at small wavevectors [p"3|. 



2.2. Kinetic energy and momentum flux densities 

The first moment and all the other odd moments of the one body density matrix are zero 
because of the reality of the wave functions. For higher even moments several definitions 
are avaible in the literature (for a review see e.g. Ref. [14]); these can be expressed again 



in terms of G(x) [|T^]. In particular, we shall concentrate on the quantity 
PJx) = (-i) 



dx 



— lim Im y^iipi \G(x)p n \ipi) 

and on its symmetrized form 



(6) 



SJx) 




p (l \x + r/2,x-r/2) 



N 



lim Im ^(ipi \p a G(x)p n - a \ ^) 



(7) 



Here we have set H — 1 and m — 1. For n = 2 these definitions give rise respectively to 
twice the local kinetic energy T(x) [Eq. (H)] and to the momentum flux density U(x) 
[Eq. (0)]. Although the semiclassical limits and the integrals of these two functions 
coincide, in the quantum limit these are different functions: the kinetic energy density 
is invoked in density-functional theory, while the momentum flux density appears in the 
equations of generalized hydrodynamics. 



The same method which was used in Sec. [2.1| to express the density profile in terms 
of the determinant of the inverse of the Green's function operator can be extended to 
deal with the operators of the form p a G{x)p n ~ a [[Ll|. Again, for low values of n one is 



reduced to evaluate with recursive methods the determinant of sparse matrices, which 
are always tridiagonal in the tails. 

The kinetic energy and the momentum flux densities obtained with this procedure 
are shown in Fig. ^| for various number of fermions in ID harmonic confinement. 
Oscillations along the profile and negative tails are found in the kinetic energy density, 
while a smooth positive profile is obtained for the momentum flux density (which 
however shows oscillations in its first derivative, see the inset of Fig. |2|). Other properties 
of these curves are that the tails of the momentum flux density profile can be expressed 
analytically in terms of the von Weiszacker surface energy (h 2 /8m)\ Vn(x)\ 2 /n(x) fTsU , 
and that an approximate form of the kinetic energy density can be found using the exact 
density profile in the local-density expression T(x) = 7c 2 h 2 [n(x)} 3 /6m fl~6| . 
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2. 3. Pair distribution function 

Higher order correlations are described in the equal-times pair distribution function. 
We focus here on how the Pauli exclusion principle is reflected in the pair distribution 
function of a gas of N noninteracting fermions in harmonic confinement. The same 
results apply to a ID Bose gas with hard-core interactions, since the boson-fermion 



mapping holds for the pair distribution function [p~3 1 . 

The longitudinal pair distribution function of a quasi- ID Fermi gas is defined as 

pi\xi, x 2 ) = n(x 1 )n(x 2 ) - F(xi,x 2 ), (8) 

where the function F(xi,x 2 ) is given by 

N 

F(x 1 ,x 2 ) = ^2 0K^i)0i(^i)0K X2 )<^( X2 ) ( 9 ) 

This function can be calculated by an extension of the Green's function method 
used in the previous sections to evaluate particle, kinetic energy and momentum flux 



densities [17|. We first rewrite Eq. (||) in terms of the Green's function G(x), 

1 - 

F(x h x 2 ) = — lim V Im(^|G(x 1 )|^>Im(^|G'(z 2 )|^) ; (10) 



we then use the property 1mA • Im_B = (l/2)[Re(A • B*) — Re(A ■ B)} to obtain the final 
expression 

F(x u x 2 ) = lim {ReTr \g n { Xi ) ■ G* N (x 2 )] - ReTr \g n ( Xi ) ■ G N (x 2 )] } , (11) 

where Gn(x) is the first N x N block of the matrix G(x). In the case of harmonic 
confinement this can be evaluated again by making use of renormalization techniques 
and recursive methods: the evaluation of F(xi,x 2 ) is reduced to the calculation of the 
determinant of pentadiagonal N x N matrices, which can be factorized into the product 
of 2 x 2 matrices. 

With this Green's function method we have calculated the longitudinal contribution 
to the equal-times pair distribution function for up to iV = 100 fermions without 
particular numerical efforts. Since this function presents in general a number of maxima 
of order iV 2 , for the sake of clarity we have shown in Fig. |3| the full result only for the 
case of iV = 4 fermions. We have plotted the pair distribution function in the center- 
of-mass R = (xi + x 2 )/2 and relative r = (x\ — x 2 ) coordinates: the effect of Pauli 
exclusion is observable in real space in the r direction as a depression of p 2 l \xi,x 2 ) at 
short distances. 



3. Dynamical properties 



One-dimensional non-interacting fermions in harmonic confinement 



6 



3.1. Collisionless regime 

The dynamic structure factor for a non-interacting Fermi gas in the collisionless regime 
is given by 

2 

/( £l -)[l - fie,)} 2ttS{u - {e j - £i )/h) . (12) 



S(k,u) = J2 [ d 3 re 



We focus here on quasi-lD excitations. The condition in this case is most stringent than 
for quasi- ID equilibrium properties (indeed the Paris experiment || satisfies the latter 
condition but not the former): the transverse excited states are not involved in the 
excitation processes when (i) the transverse component k± of the momentum transfer 
vanishes, due to orthogonality of harmonic oscillator wave functions, or (ii) the energy 
transfer u is smaller than the gap between the chemical potential and the first transverse 
excited state. 

In the quasi-lD limit Eq. (|T2| ) reduces to a one-dimensional problem and the 
transverse-state wave function factorizes out. By evaluating the overlap integral in 
Fourier space and exploiting the properties of the Hermite polynomials ||18|| , we obtain 



the expression for the dynamic structure factor as |17 

S(^h) = 2ne-^e-^ £ (| [L? (k 2 x /2)] 2 . (13) 

Here h is an integer corresponding to a single-atom excitation of h quanta of the 
harmonic oscillator, L^(x) is the i th generalized Laguerre polynomial of parameter h, 
A is the anisotropy parameter for the 3D harmonic oscillator and we have set h~ = 1, 
m = 1 and uJh Q = 1. The result for the spectrum is shown in Fig. [| at two different 
values of the transferred wave vector k x . 

A good description for the spectrum is already given by the local-density 
approximation (LDA), defined as 



S LDA (k x ,uj) = J dxn(x)S hom {k Xj uj] fi(x)) . (14) 

The spatial dependence of the chemical potential fi = k 2 F /2 = TT 2 n 2 /2 is determined 
by the relation fi(x) = fi(n(x)) and n(x) = \^2N — x 2 /ir] and the expression for 
the dynamic structure factor of the homogeneous system is Sh om {k x1 u)) = Tx/k x kp if 
1^2(^x01 < oj < uj\{k x ) and Shom(k x ,u) = otherwise, with lo\^ = k 2 /2 ± k x kF- The 
integral in Eq. (|14|) can be evaluated analytically to obtain 



2 

S LDA (k x ,u) = — 



(15) 



From Fig. [|we see that the LDA description captures the main features of the spectrum. 
The presence of the external harmonic potential significantly modifies the spectrum with 
respect that of the homogeneous system. Due to the boson-fermion mapping || the 
spectrum represented in Fig. f| is also valid for the Tonks gas. 
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3.2. Collisional regime 

We turn now to the long-wavelength limit and we investigate the analogue of sound- 
wave propagation. Due to the presence of the external confinement, the collective 
excitations are quantized. We evaluate here the spectrum and the expression for the 
density fluctuations in the linear regime. 

The equation of motion for the density profile, as can be derived exactly from the 
equations of generalized hydrodynamics |19| , takes the form 

md 2 n{x, t) = dlU(x, t) + d x [n{x, t)d x V{x, t)} . (16) 

We define the collisional regime by requiring that the momentum flux density U(x,t) 
depends locally on the particle density and the velocity field through the relation 

n(x, t) = /i 2 7r 2 n 3 (x, t)/3m + mn(x, t)v 2 (x,t)/2. (17) 

In the linear regime equations flT6|) and ( |TTD lead to a closed equation for the density 
fluctuation Sn(x,t). In the homogeneous limit we obtain a phonon dispersion relation 
with velocity c = hkp/m. 

In the case of ID harmonic confinement in order to find a discrete spectrum we 
then need to assume a generalized form of the Thomas- Fermi approximation (TFA), 
which reads 

SIl(x, t) ~ 2-^5n(x, t)=2\u- V ext (x)]5n(x, t) (18) 
on 

where t(n) is the kinetic energy density and the second equality follows from the Euler 
equation for density functional theory. We assume that Eq. (|1^) holds also outside the 
classical radius X F = ^m/i/u^: this hypothesis is necessary since the solutions of the 
linearized Eqs. (|IED and ( JT7| ) diverge at the classical boundary X F \~ 1 / 2 and thus 

we cannot obtain the dispersion relation by imposing that the solution should vanish at 
x = Xp, as in the 3D case ||20|| . By matching the solutions 8n(x)-\/\x 2 — X F \ inside and 



outside the classical radius we obtain the spectrum pi 



u = nu ho , (19) 

with corresponding solutions inside the classical radius 

= coslnareeosM) 
V 1 — z 

and outside the classical radius 

*w(z) = (ki ~ f^y i) " . (2i) 

Vz — 1 

Here we have set z = x/X F . The divergence of the solutions at x — X F is unphysical, 
and is a consequence of our approximations: the linearized TFA solutions are not valid 
around the classical turning points. 
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3.3. Non-Linear Schrddinger equation 

A possibility to go beyond the TFA is to describe the collisional Fermi gas through the 
non-linear Schrddinger equation 



ihd t $(x,t) 



^ + ^<) + fi. (lll )r 



*(*,«)• (22) 



Here, t) is normalized to the number N of particles in the trap. According to 
Kolomeisky et al. ||2"2"f , Eq. (^) describes the motion of a ID gas of impenetrable 
bosons in a long-wavelength approach. Although Eq. (|2~2"D has been shown to fail in 
the description of the dynamics of the phase of the Bose gas ||, in the linear regime 
this equation describes the collisional dynamics of the Fermi gas. Indeed, by setting 
<3>(x, t) = a/ n(x, t) exp[i(f)(x } t)] one can transform Eq. (|22|) into hydrodynamic equations 
for the density n(x, t) and the phase t) of the fluid; linearization and neglect of 
kinetic energy term —d 2 \J n(x, t)/2m^ n(x, t) leads to an expression which coincides 
with the linear form of Eqs. ([16]) and (|T7|). 

In the case of a harmonic external potential a numerical solution for the low-lying 
collective modes of Eq. (|22|) has been obtained elsewhere. The numerical results are in 
agreement with the analytical spectrum ( ]T9|) and with the main features of the density 
fluctuations profiles [Eqs. © and fl2T|)] |21[. 

It is striking that the solution of Eq. (|22| ) for the collective modes in harmonic 
confinement always yields the free harmonic-oscillator spectrum ([19]), independently of 
the strength of the interactions. In the case of the quadrupole mode, this result can be 
related to an underlying scaling symmetry of the Hamiltonian giving rise to Eq. (|22|). 
This is the one-dimensional analogue |f23[ [24]] of the hidden symmetry already noticed 



by Pitaevskii and Rosch for the 2D Gross-Pitaevskii equation | ]25[| : in the Hamiltonian 
the interaction-potential term 

H 2 7T 2 

v(x — x') = n(x)5(x — x') (23) 

3m 

scales under space dilatations x — > Xx in the same fashion as the kinetic-energy term, and 
this property in the case of external harmonic confinement leads to a closed equation 
of motion for the expectation value of the quadrupole-mode operator / = (mx 2 /2), 
expressed only in terms of the total energy E: 

3 2 t I = -Au 2 ho I + AE . (24) 

This yields uj = 2ujf lo for the frequency of the quadrupole mode. It is interesting to 



compare this value with what is obtained from the ID Gross-Pitaevskii equation |26 
namely u = ^/5/2uj ho . 



4. Summary and concluding remarks 



In summary, we have studied the equilibrium and dynamical properties of a one- 
dimensional Fermi gas under harmonic confinement. We have developed a method which 
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allows to evaluate the exact equilibrium profiles of the system for up to a large number 
of fermions. The density profiles show a shell structure which is very prominent in this 
reduced dimensionality. We have calculated the excitation spectrum of the Fermi gas 
both in the collisionless regime, where an analytic exact expression can be obtained, and 
in the collisional regime, where we have given a solution in the generalized Thomas- Fermi 
approximation. The spectra obtained in the two regimes coincide - this is the analogue 
of the coincidence of the velocities of zero and first sound in the ID homogeneous 
Fermi gas. The results here presented for the density profiles, pair distribution function 
and excitation spectrum describe as well the ID Bose gas in the Tonks regime. The 
measurement of excitation frequency of the quadrupole modes should be a sensible probe 
of the Tonks-gas phase. 
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Figure 3. Longitudinal contribution p 2 (r; R) to the pair distribution function for 
7V=4 fermions in quasi-lD harmonic confinement, in units of a^, as a function of the 
center-of-mass coordinate R and of the relative coordinate r, both in units of a^ . 
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Figure 4. Dynamic structure factor for a ID ideal Fermi gas in harmonic confinement 
(in units of w^ 1 ) as a function of fttu/E F , where E F — Nhio ho . Left panel: at k = O.lkp 
and N =500 (solid line) and 100 (bold solid line). Right panel: at k = 2k F and N =20 
(solid line) and 40 (bold solid line). Both panels show the LDA spectrum in dotted 
lines and the homogeneous result in dashes. 



